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Resume 

L'objectif de cette note est de proposer un modele mathematique simple permettant 
de comprendre I'arret de la penetration d'un flux de vapeur d'eau condensable 
sur un mur de beton, observe experiementalement (voir par exemple les differentes 
experiences d'injection de vapeur dans du beton presentees dans [11], et [7]). Un 
modele homogeneise simple d'injection dans un milieu poreux est propose, donnant 
une borne pour la position asymptotique en temps du front de penetration. 

The aim of this paper is to propose a simple mathematical model to understand 
the decision of the penetration of a stream of water vapor condensing on a concrete 
wall, observed experimentally (see for example the situations described in [11], et 
[7]). A simple homogenized model for the injection in a porous medium is proposed, 
giving a bound for the asymptotic-time position at the front of penetration. 



1 Introduction 



Un modele simple d'injection d'un fluide condensable incompressible dans un 
milieu poreux empli d'un autre fluide incompressible est presente dans [9]. 
Pour cela les auteurs ont adapte a cette situation un autre modele (cf. |8j) 
qu'ils avaient obtenu par homogeneisation a partir du systeme de Navier- 
Stokes et qui decrivait I'injection isotherme d'un fluide incompressible non 
condensable dans un milieu poreux empli d'un autre fluide incompressible. 
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L 'adaptation consistait a definir la vitesse du fluide condensable entrant comme 
somme de la vitesse qu'il aurait s'il n'etait pas condensable ( done verifiant la 
modelisation decrite dans [S]) et d'un terme de "recul" traduisant la diminu- 
tion de volume provoquee par la condensation. 

Cela leur a permis notamment de montrer que la vitesse du deplacement de 
I'interface entre les deux fluides etait decroissante et asymptotiquement nuUe 
pour les temps infinis. 

Ainsi ils purent avancer une explication pour expliquer I'arret observe experiment alement 
de I'avancee d'une vapeur condensable injectee dans du beton. Explication qui 
complete celle qui etait jusqu'alors donnee dans la litterature pour expliquer 
ce phenomene : a savoir, la formation d'un "bouchon" de condensation (voir 
par exemple ou [7]). 

Tout cela est par ailleurs developpe et detaille dans la these de Zeltz (cf. pjj) 

La demarche adoptee par ces auteurs, brievement resumee, est la suivante : 

lis considerent que la condensation provient essentiellement de la difference 
entre la pression reelle et la pression de vapeur saturante du fluide. 

Et que done cette condensation ne se fait pas uniquement au niveau de I'inter- 
face entre le fluide injecte condensable et le fluide residuel non condensable. 
Mais qu'elle se fait de fagon homogene dans toute la portion poreuse mouillee 
par le fluide injecte. 

Plus precisement, le terme utilise par ces auteurs pour modeliser ce recul est 
de la forme lineaire SV{t), V{t) etant le volume du fluide entre a I'instant t 
et S un coefficient he a la capacite de condensation du fluide considere, a la 
temperature donnee. 

L'idee precise sur laquelle repose ce choix pent s'exprimer ainsi : 

Si Pi represente la pression appliquee a la vapeur injectee et p^s la pression de 
vapeur saturante, la difference pi— p^s est constante puisqu'en condition iso- 
therme p^g est constante. Or la condensation est proportionnelle a p^— p^^ 

(cf. par exemple [S]). Et done dans ces conditions, la quantite de fluide qui se 
condense dans le milieu poreux est effectivement et uniquement proportion- 
nelle au volume du fluide deja entre. 

La question etudiee dans les paragraphes qui suivent est celle-ci : 

En gardant le point de vue qui vient d'etre decrit, le resultat obtenu dans [9] 
(deceleration progressive de I'avancee du fluide condensable injecte jusqu'a une 
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Vitesse nulle ) est-il generalisable dans le cas d'une injection non isotherme ? 
1.1 Description du milieu poreux et des autres hypotheses 



Le but principal de cette etude est d'obtenir une modelisation a la fois simple 
et fiable d'une injection de vapeur dans une enceinte en beton empli d'air 
humide. 

La structure poreuse envisagee et les hypotheses physiques ont done ete choi- 
sies en fonction. 

La modelisation de la structure poreuse est celle qui a ete utilisee dans [9j et 
dans ^ : 

Les pores ont la forme de "slits", c'est-a-dire de fentes tres fines. lis sont 
supposes paralleles et non interconnectes. La figured] schematise un tel milieu. 




I=l20mm 



Figure 1. Schematisation du milieu poreux 

Chaque slit est done un parallelipipede de dimensions e x /i x L de I'espace 
(0,5^,5^,4) 

avec e << h « L. 

L'ecoulement dans ces slits pent alors etre assimile a un ecoulement de Hele- 
Shaw, c'est-a-dire un ecoulement visqueux entre deux plans paralleles voisins. 

Ce qui permet, en procedant de la meme maniere que dans 1 10], de se ramener 
a des pores 2D (x2, X3) G Vt =]0, -L[x]0, h[ du plan (O, x^, X3). 

Un tel pore est schematise dans la figure |5] 

Le milieu poreux, et done aussi chacun des pores, est initialement rempli d'un 
premier fiuide residuel a pressions et temperatures ambiantes. Ce fiuide 
est suppose newtonien, visqueux et incompressible. Cela pent etre considere 
comme etant le cas pour de Fair en deplacement lent, car alors sa compressi- 
bilite est negligeable (pour la validite de cette approximation, on pourra voir 
par exemple P]). 
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Figure 2. Schema d'un pore. 

Le fluide injecte fi est aussi newtonien, visqueux et incompressible. II a la 
propriete supplementaire qu'il se condense lorsque sa pression depasse la pres- 
sion de vapeur saturante. C'est le cas de la vapeur d'eau. Ses pressions et 
temperatures a I'injection sont constantes et supposees tres superieures aux 
pressions et temperatures ambiantes. 

Dans la suite, la notation suivante sera systematiquement utilisee : ce qui se 
rapporte au fluide residuel fr sera indice par r et ce qui se rapporte au fluide 
injecte fi par j. 



Pour les viscosites rji et rjr des deux fluides en jeu, elles sont assez voisines 
pour pouvoir les confondre, ce qui est le cas en particulier pour Pair humide 
et la vapeur d'eau (cf. Done ces viscosites sont supposees verifler : 
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r]i{x2,xs,t) = r]r{x2,X3,t) =: r]{x2,Xs,t) dans Qx]0,T[ 



Puisque ceux-ci sont essentiellement sous phase gazeuse, les tensions superfi- 
cielles sont negligeables et il n'en est done pas tenu compte. 

De meme, les forces de gravitation peuvent etre considerees comme minimes 
pour des gaz evoluant dans un tel milieu poreux, et sont done negligees. 

Les conditions d'immiscibilite des deux fluides en jeu, de transmission a I'in- 
terface (continuites de la vitesse et de la composante normale du tenseur de 
contrainte) et d'adherence aux parois (vitesse nuUe au bord) sont celles qui 
sont habituellement retenues dans ce genre d'ecoulement de fluides newtoniens 
visqueux incompressibles. 

Elles sont identiques a celles deja utilisees dans [9]. 

Par contre, et c'est la seule difference dans les hypotheses avec celles adoptees 
dans [9] , les conditions ne sont pas isothermes et revolution de la temperature 
provoquee par la conduction de la chaleur a travers le milieu poreux est prise 
en compte. Cela sera precise par la suite. 

1.2 Developpement asymptotique et caracterisation de la vitesse 

A tout instant t > 0, le domaine fl est reunion disjointe de la partie Qr{t) en- 
core emplie du fiuide non condensable residuel et de la partie Qi(t) deja emplie 

du fiuide entrant condensable : = (t) U i7j (t) avec Q.^.{t) f] ^li{t) = ■ 

L'interface T{t) a I'instant t correspond done a ^l^(t) fl fii(t) • 

Cela est schematise dans la figure O 

La vitesse v dans Q pent etre decomposee par : 

V{X2, Z, t) = 1^{X2, Z, t) — k{x2, t) 5C2 

avec et k ainsi definis : 

- Ij^ represente la vitesse qu'on aurait pour le meme probleme mais sans 
condensation. 

- k{x2,t) est un reel positif ou nul qui correspond a la perte de vitesse pro- 
voquee par la diminution de volume qu'entraine la condensation. 

Les developpement s asymptotiques sont ecrits en fonction de e = h « L 
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Figure 3. Injection de la vapeur dans le pore. 



a partir du domaine =]0, L[x]0, e[ 



Les inconnues que constituent la vitesse du fluide, sa pression p'^ et sa 
viscosite r}'^ sont redefinies sur f2 = et dependent des variables z = 
et t, comme cela est fait classiquement dans la methode de reduction introduite 
par Ciarlet et al. (cf. [3]). 



Les developpements utilises sont precisement ceux-ci : 
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k>0 



< p''{x2,z,t) = e '^p"{x2,t) + Y,e^ '^p^{x2,z,t) 

k>0 

v^{x2,z,t) = Yj^^^^{x2-,z,t) avec v^{vl,vl) 

k>0 



A partir du systeme de Stokes et comme cela est montre dans le chapitre 4 
de [12], le systeme homogeneise obtenu en considerant uniquement les termes 
d'ordre zero est alors le suivant : 



+ =-^{*V.,*)} da,. nxlO.TI (1) 

^ + ^{V'^D + 7^{';S»} = -1^ {'l''k''(x„t)} dans nx]0,T[ 

I at oz 0x2 0x2 

avec pour conditions initiales et au bord : 



' {}*(0, z, t) = (moo, 0) dans ]0, l[x]0, T[ 
iP{L,z,t) = {uoL,0) dans ]0,l[x]0,T[ 
JqUqo dz = q sur ]0,r[ 

< 

if>{x2, z, t) = -k^{x2, t)5t2 dans ]0, L[x {0, 1} x]0, T[ 
ri%X2,z,0) = rjr sur ]0,L[x]0,l[ 
r]%0,z,t) = r]i sur ]0,l[x]0,T[ 



1.3 Comportement asymptotique de I'interface 



L'obtention du comportement asymptotique de I'interface se fait en supposant 
que la temperature est donnee par la solution de I'equation de la clialeur 
modelisant revolution de la temperature provoquee par conduction en milieu 
unidimensionnel semi infini. Nous preciserons d'abord les conditions initiales 
et au bord de la temperature et de la pression. L 'evolution de la pression de 
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vapeur saturante en fonction de la temperature est alors deduite. Ensuite est 
determinee revolution de pression soumise a la vapeur a I'interieur du milieu 
poreux. Nous definissons alors la fonction de condensation 6 : (x2; t) i — > 
S{x2;t). Cela permet de determiner la fonction k{x2,t) intervenant dans la 
definition de la vitesse. 

Finalement nous precisons la position X2{z,t) de I'interface et nous etudions 
son comportement asymptotique quand t tend vers +00. 

1.3.1 Conditions initiales et au bord 

A I'entree E = {x2 = 0}, les temperatures 6 et pressions p sont maintenues 
a partir de I'instant initial t = aux valeurs constantes 0(0, t) := 0^; et 
p(0, t) := pe- De meme, apres la sortie S = {x2 = L}, la pression est maintenue 
a la valeur constante p{L, t) : = ps- 

Enfin, a I'interieur du milieu poreux, les conditions initiales de temperature et 

de pression sont les suivantes : Pour tout X2 g]0; L] : 0(x2, 0) = 0^ et p{x2, ^) = Ps ■ 

En raison des problemes concrets qui sont modelises ici, ces valeurs verifient : 
Qe > 05- et Pe > Ps- 

1.3.2 Determination de la temperature et de sa limite asymptotique 

Les transferts de chaleur dans le milieu poreux considere sont supposes se 
faire uniquement par conduction thermique et non par convection. Cette ap- 
proximation est valide si la penetration du fiuide dans le milieu poreux est 
tres lente et si les temperatures du fiuide injecte depassent 100°C comme cela 
est clairement etabli dans I'etude parametrique resumee pages 195 et suivantes 
dans la these de Laghcha (cf. [7J). D'autres justifications de la validite de cette 
approximation se trouvent par exemple dans [10] ou encore dans [1]. 

La temperature est done supposee obeir uniquement a une equation de 
chaleur issue de la loi de Fourier. 

Cette loi est supposee de plus unidirectionnelle, puisque les flux de fluides 
dans les experiences considerees sont orientes des masses chaudes de I'entree 
vers les masses froides de la sortie. 

EUe est done du type = 011 K est une constante strictement positive 
dependant uniquement des proprietes physiques du milieu poreux considere, 
en I'occurence la densite, la conductivite et la chaleur specifique du materiau 
solide qui le constitue (cf. par exemple [6]). 
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Compte-tenu de ce qui precede, la temperature 9 est solution du probleme 
suivant : 



_ 1 ae 



9,2 - dans 0<X2; 0<t<T; 

9(0, t) =Qe> Qs pour < t < T ; (3) 
0(2^2,0) =65 pour < a;2 



Ce probleme est classique et admet pour solution fondamentale (cf. [2] page 
106) : 



et impose done la valeur suivante a Q{L,t) : 



e(L, t) = Qs + {Qe - 65) 1^1-^ /o'^ e( - x^)dxj 

Remarque 1 : Pour tout X2 € [0; L] : lim 0(x2,t) = Qe pour tout t fixe 
dans R+ : lim Q{L,t) = Qs- 

L— >+oo 

Remarque 2 : La viscosite pour un fluide donne n'etant fonction que de la 
temperature, et celle-ci etant ici fonction uniquement de X2 et de t, la viscosite 
commune des deux fluides verifie done : 

r]{x2, z,t) = r]{x2,t) dans Qx]0,T[ 



1.3.3 Evolution de la pression de vapeur saturante et de sa limite asympto- 
tique 

En supposant que la vapeur en jeu se comporte comme un gaz parfait, la 
formule de Clapeyron (cf. encore [6]) permet d'evaluer la pression de vapeur 

saturante par : = pgc^^^"'^^ 
avec : 

©0 : temperature d'ebullition de la substance a une pression po donnee, en 
degres Kelvin. 

Pvs '■ pression de vapeur saturante, dans la meme unite que po 



9 



A : une constante strictement positive dependant de la masse molaire et de la 
chaleur latente de la substance. 



En prenant pour temperature de reference Go la temperature 0^ et en notant 
lis la pression pour laquelle le fluide fi est a ebullition a la temperature Qs, 

la pression de vapeur saturante du fluide verifie done p^^ = Tr^e'^'^'^s"©'' 

oil Aj est la constante A associee a ce fluide fi. 

Compte-tenu de la Remarque 1, il vient : 

Pour tout X2 G [0; L] : lim Pvs{.X2,t) = TCse ®s ®b . 



1.3.4 Determination de la pression appliquee et de sa limite asymptotique 

La relation — 7^|??°-;7^| + t\ — = dans fixlO, T\ dans le systeme[Tl obte- 

oz oz 0x2 

nue par homogeneisation a des equations de Stokes (cf. [H] et [12] par exemple) 
permet de calculer la pression pi imposee a la vapeur (qui est dans la suite 
confondue avec sa pression essentielle po, le terme d'ordre de son developpement 
asymptotique). 

V2 verifiant le systeme [T] est donnee par (voir proposition 15 p. 52 [12]) : 



v'2ix2,z,t) = q^^^^-k%X2,t) (4) 



oil les fonctions P et R sont ainsi definies : 
^ R{x2,t) = (/o ;^o^^|;^)(/o r?o(x2^,5,t)) ~ (/o '^if^i^y 



En utilisant la Remarque 2, un calcul simple montre que : 

2 2" ^" ""V-^^) "J 2 



P^X2,z,t) = \z-\z^ et R{x2,t)-^ 



Ceci donne : 

v^{x2,z,t) = q{z-z'^)-k°{x2,t) (5) 
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3tP d dv^ 

Ainsi nous evaluons la derivee de la pression avec : = — — ^| = —2qri'^(x2, t) 

dxo dz dz 



Done la pression de la vapeur injectee verifie : 



Pi{x2,t) = -2q ri\(^^t)dC + PE (6) 

II reste a preciser la viscosite r]{x2,t) = rf{x2,t) pour preciser cette pression. 
Or elle-meme depend de la temperature suivant la loi de Sutherland (cf. [6]) : 



oil Ton note r]^^^ la viscosite du fiuide a la temperature d'entree 9^; et oil \E' 
est une constante dependant du fiuide visqueux considere. 

II en resulte que pour tout X2 G [0; L] : lim rp{x2,t) = r]^(x2,oo) = tjq 

Mais par ailleurs, la pression pi de la vapeur injectee decroit jusqu'a la valeur 
ps qui est imposee par la sortie, et reste done ensuite constante a cette valeur. 



Done, en tenant compte de cette remarque et en passant a la limite dans [6l la 
pression pi tend v€ 

definie d'apres[6] : 



pression pi tend vers une limite asymptotique Pi{x2, oo) := lim Pi{x2,t) ainsi 



Pour tout X2 G [0; L], 



(—2qriQ„X2 + Pe si X2 < ^^—^ 
" (7) 
Ps sinon. 



Remarque 3 : La pression asymptotique Pi{x2, oo) obtenue par cette modelisation 
est done essentiellement une fonction affine decroissante de la variable d'es- 
pace X2- Cela est en assez bonne correspondance avec le comportement de la 
pression observee dans les differentes experiences d'injection de vapeur dans 
du beton faites par Shekarchi et par Laghcha dans leurs theses respectives JTT], 
et Voir par exemple les figures V.29 et V.30 page 202 de [77]/ et les figures 
IV. 6 et IV 7 pages 186-187 de ^. 
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1.3.5 Determination de la fonction de condensation 5 et de sa limite asymp- 
totique 



Nous avons suppose que la condensation se faisait continuement jusqu'a I'avancee 
du front et non seulement au niveau de celui-ci. 

En effet, le phenomene de condensation de la vapeur provient essentiellement 
de I'existence d'une difference positive entre la pression imposee au fluide et 
sa pression de vapeur saturante (cf. par exemple [5). 

Et cet ecart est present sur toute la partie mouillee par le fluide condensable 
injecte, et pas seulement au niveau de I'interface avec le fluide residuel. 

Aussi le calcul de la fonction de condensation 5 est le suivant : 

Soit Pi la pression initiale de la vapeur injectee. 

La vapeur qui se condense est celle qui est necessaire pour que Pi prenne la 
valeur de la vapeur saturante Pvsi.^)- 

Plus precisement : 5{x2, t) = max (O; ^'^'''f g";/^''*^ ) 

Le passage a la limite donne pour tout X2 G [0; L] : lim 5(x2,t) := Soo{x2) 
oil 

[ ^ _ ^ge^''^"^^ ^ PE-PS 

Soo{x2) = \ P-;^T-7, '-2..e, (8) 

^ sinon. 

I PS 

Remarque 4 : Comme cette etude s'interesse avant tout au comportement 
asymptotique lorsque le temps t tend vers I'infini, t est suppose dans la suite 
assez grand pour pouvoir approximer 6{x2,t) par Sooix2). 

1.3.6 Determination de la vitesse de recul de I'interface 

Pour un z fixe dans [0; 1] et pour un instant t fixe, soit X2 = X2{z, t) la position 
correspondante de I'interface entre le fiuide residuel et le fiuide condensable. 

Suite a la condensation du fiuide condensable present sur I'horizontale de cote 
z entre X2 = et 2:2 = X2{z, t) , cette interface subit une reculee dont la vitesse 

est ko{z,t) = Jq^^^'^^ S{x2,t)dx2 
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1.3.7 Comportement asymptotique de X2{t,z) 

Comme dans le cas isotherme envisage dans |9] et [12], montrons que la posi- 
tion de I'interface est bornee pour tout temps. 

D'apresOla vitesse est donnee par : 

fo(x2, z, t) = wo{z) - ko{z, t), 

la position X2{z, t) a I'instant t de I'interface entre les deux fluides qui etait en 
M(x2 = 0, z) a I'instant t = verifie done : 

X2{z, t) = Wo{z)t - Jq ko{z, T)dT = Wo{z)t - /o /o''"^'''*'' 6{X2, t)dX2dT 
Et puisque pour les grands temps nous avons : 

5{x2,t) ^ 5^{X2) > 1 - ^ > 

Ps 

la position X2{z,t) de I'interface est majoree par la solution de I'equation 



6 = Wo{z)t - 1 - / C2dT 

Ps / -^0 



c'est-a-dire par : 

wo{z)t 



PS 



Et done dans ce cas nous avons : 

Wo{z)t 



lim X2{z, t) < lim 



1+ 1 - t 



soit finalement : 



lim X2{z, t) < "^f. , , . < +00. (9) 



PS 



1.3.8 Conclusion 
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En se plagant dans des conditions non isothermes, avec un modele tres simple 
pour la temperature, le comportement global deja obtenu dans des conditions 
isothermes dans (cf. [9j et [12\) est maintenu : 

S'il y a condensation du fiuide entrant, la vitesse de I'interface entre le fiuide 
entrant et le fiuide residuel decroit progressivement vers et la position du 
front tend vers une position asymptotique finie quand le temps t tend vers 
I'infini. 
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